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^ : Abstract 

We study the Josephson effect through a one-dimensional system of inter- 
tH- ■ acting electrons, connected to two superconductors by tunnel junctions. The 

CN , interactions are treated in the framework of the one-channel Luttinger model. 

1^^ I At zero temperature, the Josephson critical current is found to decay alge- 

^D ' braically with increasing distance between the junctions. The exponent is 



a 



proportional to the strength of the Coulomb interaction. If the Luttinger liq- 
uid has a finite size, the Josephson current depends on the total number of 
electrons modulo 4. These parity effects are studied for the ring, coupled ca- 
pacitively to a gate- voltage and threaded by a magnetic flux. The Josephson 
'^ I current changes continuously as a function of the gate voltage and stepwise as 

bi • a function of the magnetic flux. The electron-electron interaction introduces 

CJ ■ qualitatively new features compared to the non-interacting case. 
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Recent technological developments in the fabrication of semiconductor- superconductor 
(S-Sc) interfaces made the observation of supercurrent in S-Sc-S junctions possiblelll. This 
progress was accompanied by a number of appealing theoretical predictions. The Josephson 
current through a narrow constriction in a high-mobility two-dimensional non-interacting 
electron gas (2DEG) should be quantizedB. Interference effects of two electrons entering 
the normal region from a superconductor are important, especially in mesoscopic samplesQ. 
This interference influences the Josephson current through the normal region. 

It is well-known that electron-electron interactions affect the transport properties of 
mesoscopic devices as well. The interaction manfests itself in the single charge effectsj. 
In particular, the Josephson curent through SSS and SNS systems can be modulated by a 
gate-voltageB. This has been described in terms of a phenomenological capacitance model. 
However, in low-dimensional semiconductor nanostructures, with a low electron density, 
electron-electron interactions should be treated microscopically. As a result, the transport 
properties of a (quasi-) one-dimensional quantum wire show deviations from Fermi-liquid 
behaviouiO. 

In the present paper we analyze the Josephson current through a Luttinger liquid (LL). 
Specifically, we consider two geometries which can be realized experimentally: A long wire 
with tunnel contacts to two superconductors at a distance d (Fig. |I|a), and a ring with 
circumference L and tunnel contacts at a distance L/2 (Fig. |l]b). Various aspects of trans- 
port in mesoscopic systems (parity effects and interference combined with electron-electron 
interactions) and their interplay can be studied in this system. The aim of this paper is 
twofold. First, we want to study the influence of the Coulomb interaction on the Josephson 
critical current. Second, we would like to see how the parity effects which are present in the 
ringpil3 manifest themselves in the dependence of the Josephson current on magnetic flux 
and gate-voltage. The coupling of bulk superconductors to a one-dimensional (ID) system 
of interacting electrons was studied by FisherQ. In contrast to our work, a chiral LL of 
spin-polarized electrons was considered and tunneling with spin flip was analyzed. 

The systems under consideration (Fig. ^ can be described by the Hamiltonian H = 
Hsi+Hs2+Hl+Ht- Here, Hsi, Hs2 are the BCS-Hamiltonians for the bulk superconductors 
with gap Ai, A2 respectively, kept at a phase difference X = Xi ~ X2- The (ID) electron 
system is described by the HamiltonianQ 









(1) 



It is written as a sum of the contributions from the spin (j = a) and charge (j = p) degrees 
of freedom. We use the standard notations for the parameters gj of the interaction strengths 
and the velocities Vj of spin and charge excitations. We also introduce the bosonic fields 
'^s = 4>p^ -500- and 6s = 6p + s9a for spin up (s = 1) and down (s = —1) fermions. These 
fields obey the commutation relation [(/>s(a;), 9s'{x')] = (iii/^) siga{x' — x)6s,s'- The tunneling 
is assumed to occur through two tunnel junctions at the points x = and x = d, 

Ht = Y. Ti^i,s{x = 0)^L,.(x = 0) + T2^^,(x = d)^LA^ = d) + (h.c). (2) 

s 

The tunnel matrix elements Ti 2 can be related to the tunnel conductances G'i_2 of the 
junctions: Gi = {4iTe^/h)NL{0)N^{0)T^, where Nl{0) = l/Tchvp, Ni is the normal density 
of states of a superconductor {i = 1, 2), and vp is the Fermi velocity. 



The fermionic field operators ^ can be expressed in terms of spin and cliarge degrees of 
freedom B: 



^L(^,^) 



m,odd 



exp {imkpx} exp {imOs} exp {i4>s} 



(3) 



where kp is the Fermi-wave vector and po,s = Ng/L = kpJi2Tx) is the average electron 
density for one spin direction. If the LL is confined to a ringll3, threaded by an Aharonov- 
Bohm flux $, the twisted boundary condition ^!\ ^{x + L) = exp {— 27rz$/$o}^L s(3^) can 
be imposedli^. Here, $o is the normal flux quantum hje. The fields and are decomposed 
in terms of bosonic fields and topological excitations^ in the following fashion: 



ej{x) 



,-(x)+^° + 7rMj(x/2L), 



ct>,{x) = 0,(x) + cf + 7r(J, - A5j,p^/%){x/2L). 
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Here, 9j and 0j are the non-zero modes 



e,{x) 



0,-(x) 




sign(g)e^^^'(^],, 



h,-q)^ 



^iqx 
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where hj^q, bj^ are Bose operators. Mj and Jj denote the topological excitations. They are 
related to the usual topological excitations for fermions with spin s: Mg = (1/2) [Mp + 
sMa-] and Jg = (1/2) [Jp + sJo-]- Using the topological constraints for Mg and Jg given 
in Ref.a, we find constraints for Mj and Jj: (i) The topological numbers Mj and Jj are 
either simultaneously even or simultaneously odd; (ii) when Ng is odd the sum Mp ± Mo- + 
Jp ± Jo- takes values ..., —4,0,4, ..., when A^^ is even the sum Mp ± M^- + Jp ± J^ takes 
values ..., —6, —2, 2, 6, .... Here, the number of electrons A^^ (A^s=i= A^5=-i) determines the 
linearization point kp = nNg/L of the original electron spectrumM. 

The Hamiltonian is decoupled in the non-zero modes and the topological excitations: 



HL = hY: 



Y.^1 



TTVi 



\m,Aj + -jT 



J=P.o" 197^0 



AL 
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4^i,p/<i 



+ -{Mj 
9j 



"^^pU 



(6) 



We introduced the flux frustration /$ = $/$0) as well as the parameter fp = {gpL/4TTVp)Afi, 
which is related to the difference A/i between the electro- chemical potential fi of the super- 
conducting electrodes and the Fermi energy Ep^ of the ring. For non-interacting electrons 
the latter is just the energy {hkpY /2M in the linearization point. Generally, the reference 
point A/i = is defined from the requirements that for $ = there are 2A^g electrons in the 
ground state (A^i = A^-i) and the energies to add and remove the electron(s) to/from the 
system are the same. The electro-chemical potential /i can be controlled by a gate-voltage. 
The stationary Josephson effect can be obtained by evaluating the phase- dependent 
part of the free energy T{x)'- The Josephson current is given by Ij = —{2e/h)dJ-'/dx- We 
expand JF = — (1//5) InZ where Z = Tt exp{—j3H} in powers of Hp] the lowest order con- 
tribution arises in 4*'^ order. It is represented by the diagram in the inset of Fig. ^, which 



has a clear physical meaning: the Josephson effect occurs through a process which trans- 
fers a Cooper pair from superconductor 5*2, described by the anomalous Green's function 
Fs2{t3,T4) through the second junction with an amplitude (Tg*)^. Subsequently, the two 
electrons propagate as a Cooperon through the LL from contact 2 to 1. This propagation 
is determined by the Cooperon propagator 11(0, d; ri, ..., r4). Finally, both electrons tunnel 
through the first junction (amplitude Tf) and enter the superconductor SI as a Cooper 
pair, characterized by the anomalous function Fgi{Ti,T2). If the distance d between the 
junctions is much larger than the coherence length in the superconductor, the characteristic 
energies Tivp/d of the electrons propagating through the ID system are much less than the 
energies A of excitations in superconductors. At low temperatures (fc^T ^ A), a generic 
process consists of fast tunneling of two electrons from the superconductor into the ID sys- 
tem (|ri — r2| ~ ^/A) and their slow propagation through the LL (|ri — TsI ~ d/vp)- The 
phase-dependent part of the free energy then simplifies to 



T{x) = -27r2iVi(0)iV2(0)3fJe[T2(T2*)V 



■JX 



dTli{d]T)], 



(7) 



whith the Cooperon 



n(c?,r) = (T.^^,+ (0,0)^i,_(0,0)^i,_(t/,r)^l , (rf,r) 



(8) 



where the average is taken over the eigenstates of Hl- The evaluation of (H) with the help 
of bosonized field operators like (|^) is straightforward. 

Wire geometry. For an infinitely long wire (Fig. [l|a), the topological excitations play no 
role since their energies ~ t^Vj/L are vanishingly small. The Cooperon propagator n((i, r) 
is given by 



n^((i, r) 



pI 



y^ ^i{ni+n2)(kpd+ri) ^ 
ni,n2,odd 



a 



d^ + vIt^ 



go-(™i-n2)^/16 



a 



d^ + f ^r^ 



l/gp+9p{ni+n2f/l& 



(9) 



where e*'' = [{d+iv pT) / {d—iv pT)]^/"^ , and a is a cut-off parameteiiB of the order l/kp- If ni-|- 
n2 7^ 0, the rapid oscillations related to e*("i+"2)(feFrf) j^iake the Josephson term vanishingly 
small if the tunnel junctions are large compared to the Fermi wavelength. For ni + ^2 = 
the leading contributions to the correlation function correspond to ni — n2 = ±2. For the 
spin-independent electron-electron interaction we should fix gfo- = 2 and Vp = {2/ gp)vp. For 
the Josephson current at zero temperature we find 



2Tievp G1G2 

d {Aeyhy 



;F^{gp,d) sm{x), 



(10) 



kpd 



2/9P-1 rdx 



71" vT+x^ 



{Fw = 1 for noninteracting electrons). 

Ring geometry. For the ring (Fig. |l]b) we get 



1 + (2x/^,)2 



1/9P 
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Here, (...)j,a/ means evaluation with respect to the ground-state configuration of the topo- 
logical excitations Jj, Mj of the LL subject to the topological constraints. 

In Fig. ^ the dependence of the critical current on Qp for the two geometries is shown. 
The Josephson current is suppressed by the repulsive interaction, Ij oc d~'^/^'' . We estimate 
the magnitude of the current, using typical numbers from the experiment of Mailly et al.lij, 
to be of the order of several nA in the non-interacting case. 

We discuss now the flux and gate-voltage dependence in the ring geometry. The Joseph- 



son current Eq. (|T3|) depends on the gate- voltage via the parameter /^. The flux does not 
enter expicitly into Eq. (jT^). However, the Josephson current depends on /$ (and on /^) 
implicitly via the topological numbers {Jp, J^, Mp, M„). Therefore, one can expect that the 
Josephson current changes stepwise as a function of the flux (the jumps correspond to the 
change in the topological numbers) and shows a continuous dependence (with jumps) on the 
gate- voltage (see Figs. |], ||). 

Without loss of generality we restrict the further consideration by odd values of A^^"- We 
start from the non-interacting case gp = 2. The ground-state conflgurations {Jp, J„, Mp, M„) 
are displayed in Fig. ^. Each state in the ring below the chemical potential /i of the 
superconductors is occupied by two electrons with opposite spin. Hence, the topological 
numbers M„ and Jo- are always zero. The increase of the flux shifts down (up) the energy 
levels of the electrons moving in the positive (negative) direction along the ring. Increasing 
fp correponds to a uniform shift of all energy levels down with respect to the chemical 
potential /i. A pair of electrons tunnels into (out of) the ring each time when an empty 
(fllled) energy level crosses the chemical potential. This leads to a change in the topological 
number Mp, AMp = 2(— 2). Simultaneously, the topological number Jp changes, AJp = ±2, 
since the system acquires orbital momentum. With the increase of the flux $ by one flux 
quantum, a pair of electrons enters the positive branch of the spectrum and (for a different 
value of $) another pair leaves the negative branch. For each process AJp = 2. Accordingly, 
the Josephson critical current Ij^c = niax^|/j(x)| (Fig. ^ changes stepwise with the flux 
making two jumps per period $o- One of the jumps occurs at /$ = — /^ + 1/2 in Fig. |. 

For fp, = 0, the chemical potential of the superconductors /i is in the middle of the gap 
between the last occupied and flrst unoccupied energy levels {Ep^Q = jj). For this reason, 
the two processes mentioned above occur at the same point /o- = 1/2. One can say that 
a pair of electrons jumps from the negative branch to positive branch (e.g. the transition 
(0, 0, 0, 0) — *> (4, 0, 0, 0) in Fig. ^a) and the number of electrons in the system is conserved. 

The situation changes drastically for interacting electrons. The interaction lifts the 
degeneracy of the single particle states with different spin and the system acquires additional 
stiffness with respect to the change in the electron number (for the repulsive interactions gp < 
2 the coefficient 2/gp of the term with the particle number Mp is larger than the coefficients 
of the other topological terms in Eq. (BI)). For this reason, the electrons tunnel one by one 



into the ring and new ground-state configurations (e.g. {Jp, J„, Mp,M„) = (1,1,1,1) and 
(3,-1,11)) with an odd number of electrons arise (see Fig. ^a). In particular, the parity 
of the electron number on the ring can be changed not only by changing the gate-voltage, 
but also by changing the flux. Depending on the value of the gate-voltage the increase of 
the flux by one flux quantum is accompanied by the tunneling of two, one or zero electrons 
into the ring and the same number of tunneling events out of it. Accordingly, the Josephson 
current (Fig. ^) shows four, two or zero jumps per period $o- 

The size of the intervals of the gate- voltage, where the electron number (characterized by 
Mp) does not depend on the flux, increases with the increase of the repulsive interaction. In 
the limit of strong interaction the number of particles is fixed for any value of the gate- voltage 
except for small intervals near the points /^ = 1/8 + n/4 with integer mB. 

We have also found that for some ground-state configurations, the Josephson current 
changes sign and the ring acts as a vr-junction. This happens, for example, in the regions 
with Jp = 2 mod(4) both in the interacting and in the non-interacting cases (e.g. in the 
region (2,0,2,0) in Figs. ^, ||a). Such behaviour can be tested through the interference 
pattern of a SQUID, consisting of the ring and a standard Josephson junction. A more 
detailed analysis of all these features will be presented elsewherellll. 

In this letter we studied the Josephson effect through a ID system of interacting electrons: 
a quantum wire and a ring. The Josephson current was found to be supressed by the Coulomb 
interaction. The interaction qualitatively modifies the parity effects in the ring and effects 
the dependences of the critical current on the flux and gate- voltage. An anomalous sign of 
the Josephson current was found in some range of parameters. 
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FIGURES 

FIG. 1. The geometries discussed in the text: (a) one-dimensional wire connected to two 
superconductors by tunnel junctions, separated by a distance d. (b) Ring with circumference 
L, threaded by a magnetic flux $. It is connected to two superconductors by tunnel junctions, 
separated by a distance L/2. 



FIG. 2. Critical current as a function of gp for the ring (solid line, normal- 

r(0) 

r(0) 



ized to /jj = (47refF/L)(GiG2/(4e^/?i)^)) and the wire (dashed line, normalized to 



Ij^ = {irevF / d){GiG2/ {4:6 /h) )) . Inset: Lowest-order contribution to the Josephson effect. 

FIG. 3. (a) Ground-state configurations for topological quantum numbers {Jp, J^, Mp, Ma-) as 
a function of gate- voltage (/^) and magnetic flux (/$), for gp = 2. (b) Critical current at T = 
for (7p = 2 as a function of /$ and fp. 

FIG. 4. The same as Fig. 3, but for gp = 1.75. 



